202208.00052v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


tów 


ARID LAND GEOGRAPHY 


融合 影像 在 土壤 盐 江 化 监测 中 的 
适用 性 人 研究 
一 一 以 渭 干 河 一 库 车 河 三 角 洲 绿 》 


RIF"? Tm”, IFLA", 金 晓 叶 HR Iz -F 4 
(1. 新 疆 大 学 地 理 与 遥感 科学 学 院 , 智 慧 城市 与 环境 建 模 自治 区 普通 高 校 重点 实验 室 ,新 疆 乌鲁木齐 ”830046; 
2. 新 疆 大 学 绿洲 生态 教育 部 重点 实验 室 ,新 疆 乌鲁木齐 ”830046; 3. 自然 资源 部 中 亚 地 理 信息 技术 开发 应 用 工 
程 技术 创新 中 心 ,新 疆 乌鲁木齐 ”830046; 4. 新 疆 阿 克 苏 地 区 渭 干 河流 域 管理 处 ,新 疆 阿克苏 ”842000) 


MODIS 和 Landsat 时 空 


摘 要 : 土壤 盐 渍 化 的 遥感 监测 依赖 于 高 时 空 分 辨 率 影像 ,但 受 经 费 预 算 、 卫 星 回访 周期 及 天 气 的 
影响 ,高 时 空 分 辩 率 的 遥感 影像 较 难 获取 ,这 就 限制 了 根据 采样 时 间 来 获取 对 应 时 期 遥感 影像 进 
行 土壤 盐 渍 化 监测 反 演 的 应 用 。 为 此 ,提出 融合 MODIS 和 Landsat 影像 生成 高 时 空 分 辩 率 影像 来 
提取 土壤 盐 渍 化 信息 ,为 时 空 影像 进行 土壤 盐 渍 化 监测 研究 提供 数据 参考 。 以 渭 干 河 一 库 车 河 绿 
洲 为 研究 区 ,利用 增强 型 时 空 自 适应 融合 率 反 射 模型 (Enhanced spatial and temporal adaptive reflec- 
tance fusion model,ESTARFM) 和 灵活 的 时 空 融合 模型 (Flexible spatiotemporal data fusion ,FSDAF ) ， 
对 MODIS 和 Landsat 影 像 进行 时 空 融合 ,并 基于 融合 影像 数据 构建 了 关于 土壤 电导 率 (EC) 的 随机 
森林 (REF) 预 测 模型 ,对比 分析 时 空 融合 影像 应 用 于 土壤 盐 渍 化 遥感 监测 的 适用 性 。 结 果 表 明 :ES- 
TARFM 融合 影像 的 特征 波段 反射 率 与 Landsat 验证 影像 对 应 波段 反射 率 一 致 性 决定 系数 忆 (Red)= 
0.8066 R° (SWIR2)=0.8444; FSDAF 融合 影像 的 特征 波段 与 Landsat 验证 影像 对 应 波段 反射 率 一 至 
a =0.6999 R° (SWIR2)=0.7493 ; 基于 ESTARFM 融合 影像 构建 的 EC 值 预测 模型 
精度 最 高 ,R=0.9268, 基于 FSDAF 融合 影像 构建 的 EC 值 预 测 模 型 精度 良好 ,R=0.8987, 基于 验证 
影像 构建 的 EC {H Hi il 4% A R°=0.9103; ESTARFM 模 型 的 融合 精度 高 于 FSDAF 模 型, 并且 基于 融合 
影像 构建 的 EC 值 预测 模型 效果 良好 。 
X 键 词 : 时 空 融合 ; 随机 森林 ; 土壤 电导 率 ; TRAR 


文章 编号 : 

土壤 盐 涡 化 一 般 指 易 溶 盐分 在 土壤 中 含量 不 硬件 技术 条 件 等 的 限制 ,高 时 空 分 辨 京 的 各 感 影像 
断 积累 的 现象 或 过 程 。 发 生 盐 渍 化 的 耕地 土壤 不 依旧 很 难 获 取 , 当 前 常用 的 遥感 影像 数据 仍 存在 
利于 作物 生长 ,使 农作物 减产 或 绝收 ,对 粮食 安全 “WFE” PA, Landsat 等 遥感 卫星 影像 虽然 
造成 危害 ,同时 会 带 来 水 土 流失 土壤 侵蚀 Ew ”空间 分 辩 率 较 高 ,但 其 时 间 分 辩 率 较 低 , 且 光 学 般 
漠 化 等 一 系列 的 环境 恶化 问题 ,给 生态 环境 以 及 物 感 影像 易 受 到 云层 覆盖 等 大 气 条 件 的 影响 。MODIS 
种 多 样 性 带 来 威胁 。 定 量 有 效 预 测 土壤 盐 渍 化 情 a 周期 短 ,可 获取 高 时 间 分 辩 率 的 
况 ,掌握 土壤 盐分 分 布 规律 ,对 于 当地 上 土壤 盐 渍 化 ”遥感 影像 ,但 其 空间 分 辩 率 较 低 。 为 此 ,提出 多 种 


的 治理 具有 重要 意义 。 
随 着 空间 信息 技术 的 迅速 发 展 ,遥感 技术 在 土 
壤 盐 分 预测 中 应 用 前 景 广阔 。 然 而 , 受 预 算 经 费 和 
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多 源 遥 感 数据 融合 技术 ,以 获取 兼 具 高 时 间 和 高 空 
间 分 辨 率 的 遥感 影像 ,来 满足 大 范围 .高 精度 ,快速 
变化 的 地 表 监 测 与 大 气 环境 遥感 的 应 用 需求 。 目 
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前 较为 流行 的 多 源 遥 感 数据 融合 算法 包括 :基于 时 
空 变 化 滤波 的 融合 方法 .基于 空间 信息 分 解 的 融合 
方法 .基于 机 器 学 习 的 融合 方法 和 组 合 型 的 融合 方 
法 ”。 其 中 基于 权重 滤波 的 时 空 融合 模型 原理 简 
单 ,发 展 成 熟 , 骨 型 算法 包括 时 空 自 适应 的 反射 率 
融合 模型 (Spatial and temporal adaptive reflectance 
fusion model, STARFM ) .增强 型 时 空 自 适应 融合 率 
反射 模型 (Enhanced spatial and temporal adaptive re- 
flectance fusion model , ESTARFM ) 等 "1 ,ESTARFM 
作为 改进 的 时 空 自 适应 反射 率 融 合 模型 ,提升 了 该 
类 方法 对 于 异 质 性 地 表 覆 盖 区 域 影像 的 融合 效果 
以 及 对 不 同时 相 影 像 地 物 类 型 变化 的 预测 能 力 “。 
而 组 合 型 的 时 空 融 合 方法 由 于 综合 了 其 他 不 同类 
型 融合 方法 的 优势 , 较 全 面 地 应 对 了 各 类 时 空 遥 感 
影像 融合 方法 的 不 足 , 通 常 具 有 较 好 的 融合 效果 ， 
因而 具有 更 广泛 的 应 用 前 景 。 灵 活 的 时 空 融合 模 


监测 中 的 应 用 提供 科学 支持 和 参考 。 
1 数据 与 方法 


1.1 研究 区 概况 

渭 干 河 一 库 车 河 三 角 洲 绿洲 "位 于 新 疆 南 部 
的 塔里木 盆地 中 北部 ,隶属 阿克苏 地 区 ,辖区 范围 
包括 沙 雅 县 . 库 车 县 .新 和 县 。 根 据 实 地 采样 区 域 确 
定 人 研究 区 边界 地 理 坐 标 为 81°37’~83°59'E,41°06'~ 
42°40'N, 海 拔 892~1100 m, 总 面积 523.76x10? km’, 
研究 区 属于 典型 的 大 陆 性 暖 温带 干旱 气候 ,蒸发 强 
Bl ,降水 分 布 不 均 ,山区 多 年 平均 降水 量 为 243.0 
mm 平原 区 多 年 平均 降水 量 为 46.5 mm ,平原 区 年 
平均 蒸发 量 为 1374 mm。 渭 干 河 一 库 车 河 三 角 洲 绿 
洲 土 地 利用 类 型 主要 包括 农田 、 草 地 、 林 地、 芒 地 、 
盐 涡 地 等 ;主要 土壤 类 型 为 潮 土 . 草 甸 土 ,此 外 还 分 
MAREE .沼泽 土 和 盐 土 等 多 种 类 型 土壤 。 在 平 


型 (Flexible spatiotemporal data fusion, FSDAF) ™ 组 
合 了 基于 空间 信息 分 解 与 时 空 变 化 滤波 的 2 类 融合 
方法 ,并 融入 了 空间 插值 技术 ,提高 了 上 述 算法 在 
异 质 性 地 表 歼 盖 区 的 融合 效果 ,并 提高 了 其 对 地 物 
类 型 变化 的 预测 能 

近年 来 ,国内 外 学 者 基于 遥感 影像 反 演 土壤 盐 
分 的 有 效 方法 ,主要 是 利用 多 光谱 影像 和 经 典 统 计 
分 析 进 行 研 究 建 模 ”。 而 土壤 作为 具有 高 度 变异 的 
时 空 连续 体 ,土壤 发 育 过 程 复 杂 ,部 分 学 者 在 建 模 
过 程 中 采用 Pearson 相关 分 析 进 行 特征 变量 优选 , 忽 
略 了 环境 因子 对 土壤 的 非 线 性 作用 ”, 这 对 土壤 盐 
渍 化 反 演 精度 产生 影响 。 随 着 数据 挖掘 和 机 顺 学 
习 技 术 迅 速 发 展 , 它 能 够 有 效 解决 土壤 与 环境 因子 
之 间 非 线性 问题 ,在 土壤 盐 渍 化 反 演 中 具有 深远 的 
应 用 价值 。 

时 空 融合 算法 目前 已 被 广泛 应 用 于 生态 水 
SC) RB IC ,城市 温度 ”洪水 监测 ”等 遥感 
应 用 研究 ,但 是 在 土壤 盐 渍 化 监测 中 的 应 用 研究 仍 
十 分 欠缺 。 本 文 以 渭 干 河 一 库 车 河 三 角 洲 绿洲 为 
研究 区 ,采用 ESTARFM 模型 和 FSDAF 模 型 对 MO- 
DIS 和 Landsat 影像 进行 时 空 融合 ,根据 实际 采样 时 
间 ,获取 与 采样 时 间 相 对 应 的 高 时 空 分 辩 影 像 数 
据 , 并 基于 时 空 融合 影像 数据 构建 关于 土壤 电导 率 
(EC) 的 随机 森林 (RF) 预 测 模型 ,探究 分 析 时 空 融合 
影像 应 用 于 土壤 盐 渍 化 遥感 监测 的 建 模 效 果 和 反 
演 精 度 ,为 时 空 融合 影像 在 干旱 区 土壤 盐 渍 化 遥感 


原 区 中 下 部 地 势 较 平坦 ,地 下 水 位 较 高 ,土壤 下 层 
构成 物 颗 粒 细 ,透水 性 差 ,蒸发 作用 强烈 ,导致 盐分 
随 水 运动 ,积累 于 地 表 造 成 土壤 盐 渍 化 ,灌区 内 盐 
渍 化 面积 已 达 50% 以 上 ,其 中 严重 盐 渍 化 面积 达 
30% , 盐 渍 化 土壤 pH 为 8 左右 , 含 盐 主 要 成 份 为 所 
化 物 。 因 此 ,选择 地 区 作为 研究 区 具有 较 好 的 代表 
性 ,对 改善 生态 环境 及 农业 生产 的 发 展 具 有 重要 意 
X ,人 研究 区 采样 点 分 布 如 图 1 所 示 。 
1.2 数据 来 源 与 处 理 

实测 数据 :2018 年 7 月 11 一 7 月 23 日 对 研究 区 
开展 野外 调查 ,用 五 点 法 采集 土 样 ,在 各 采样 点 30 


口 研究 区 边界 


图 1 研究 区 采样 点 分 布 图 


Fig. 1 Distribution of sampling points in the study area 
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mX30 m 的 样 方 内 使 用 木 铲 采集 $ 个 土壤 样品 ,采样 
深度 为 10 em ,并 在 现场 进行 混合 以 创建 代表 性 的 
复合 土壤 样本 ,使 用 误差 小 于 5 m 便 携 式 GPS(Uni 
Strong G120) 记 录 每 个 采样 点 的 位 置 。 待 土壤 样品 
自然 风干 后 去 除 杂 质 ,过 2 mm (10 H FL , PKA 20 
g 土 壤 样 品 与 100 mL 去 离子 水 配置 成 水 土 比 5:1 的 
土壤 悬 间 液 ,震荡 30 min 后 静 置 24 ,在 室温 25 C 
下 提取 渗 滤 液 ,采用 配备 有 复合 电极 (Tetra Con 
925) 的 数字 多 参数 测量 设备 (Multi 3420 Set B, 
WTW Gmb H, 德 国 ) 测 量 EC 值 ,最 终 获得 有 效 样本 


数量 为 69 个 。 
遥感 影像 :Landsat 影 像 来 源 于 谷歌 地 球 云 计算 


(Google Earth Engine , GEE) 平 台 的 Tl_SR 数据 集 ， 
MODIS 影像 来 源 于 GEE 平 台 的 MOD09A1 (全 球 
500 m 地表 反 射 率 8 d 合成) 地 表 反 射 率 数据 集 
(L38d)。 遥 感 影像 数据 信息 见 表 1。 

数据 已 进行 辐射 矫正 和 大 气 矫正 ,结合 实地 采 
样 时 间 和 云 量 (<10%) ,选用 2018 年 6 月 的 一 对 MO- 
DIS, Landsat 影像 与 2018 年 8 月 的 一 对 MODIS、 
Landsat 影像 进行 时 空 融 合 , 选 取 2018 年 7 月 20 H 
的 Landsat 影 像 作为 验证 影像 ,作为 融合 结果 影像 的 
比 对 参考 。 所 有 影像 数据 包括 红 光 (Red) .蓝光 
(Blue) 、 绿 光 (Green) 、 近 红外 (NIR)、 中 红外 1 
(SWIR1) .中 红外 2(SWIR2)6 个 特征 波段 ,并 在 谷 
歌 云 计算 平台 对 数据 进行 合成 、 去 云 .裁剪 的 预 处 
H, AJH ENVI 5.3 软件 对 MODIS 影 像 进行 重 投影 ， 
并 重 采样 为 空间 分 辩 率 30 m 的 数据 ,与 Landsat 影 
像 分 辨 率 保持 一 致 。 
1.3 研究 方法 
13.1 EAHA ESTARFM 模型 是 在 STARFM 
算法 的 基础 上 进行 改进 的 一 种 经 典 时 空 融合 模型 ， 
ESTARFM 模型 融合 质量 高 ,具有 较 广 泛 的 应 用 。 
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本 文通 过 ENVI 5.3+IDL 8.5 实现 模型 ESTARFM 的 
时 空 融 合算 法 过 程 ,主要 步骤 包括 :基于 每 个 参考 
时 刻 内 的 高 空间 分 辩 率 图 像 (Landsat 影像 ) 选 取 相 
似 像 元 ;相似 像 元 权重 计算 ; 预测 时 刻 了 ,。 目标 像 元 
值 计算 。 其 中 ,中心 像 元 重建 的 表达 式 为 : 
R apa D) =T R Eao T) +t 
T R(X gaans Ta) 
式 中 :Ri 为 预测 时 相 的 高 分 辩 率 融合 数据 ; 
Goyom) 为 预测 时 刻 ( 7 ) 的 中 心 像 元 值 ; R, ~ R, 
为 2 个 基础 时 相 的 高 空间 分 辩 率 数据 ; T, a T, 为 2 
个 基础 时 相 数 据 的 时 间 权 重 ; w 为 滑动 窗口 的 大 
小 ,本 实验 中 模型 参数 滑动 窗口 大 小 w 为 20， 类 别 
数 class 设 置 为 4。 
本 文通 过 ENVI 5.3+IDL 8.5 实现 模型 FSDAF 的 
时 空 融合 算法 过 程 ,在 影像 基 对 数 上 与 ESTARFM 算 
法 保持 一 致 ,选用 一 对 影像 进行 时 空 融 合 。FSDAF™" 
综合 了 混合 像 元 分 解 和 加 权 函 数 对 方法 ,能够 比较 
好 的 对 地 物 类 型 变化 的 区 域 进 行 预测 。 该 算法 利用 
初始 时 刻 的 已 知 高 . 低 分 辨 率 影像 预测 时 刻 i 低 
分 影像 ,融合 得 到 立时 刻 高 分 影像 。FSDAFS 融合 模 
型 主要 包括 6 个 步骤:(1) 进行 ,时刻 的 高 空间 分 辩 
率 影 像 (Landsat 影像 ) 分 类 ;(2) 估计 从 时刻 到 
时 刻 高 时 间 分 辩 率 影像 (MODIS 影 像 ) 每 个 类 别 的 时 
间 变 化 ; (3) 计算 六 时 刻 的 Landsat 影 像 使 用 类 别 尺 
度 的 时 间 变 化 和 计算 MODIS 影像 每 个 像素 预测 的 
残 差 ; (4) 利用 i 时 刻 的 MODIS 影像 ,使 用 薄板 样 
条 函数 预测 高 空间 分 状 率 影像 (Landsat 影像 ); (5) 
基于 薄板 样 条 函数 的 残 差 分 布 计 算 ;(6) 利用 邻 域 
言 息 得 到 Landsat 影 像 的 最 终 预 测 。 
1.3.2 预测 模型 构建 ”REF05 是 一 种 统计 学 习 理 论 ， 
它 使 用 自 举 重 采 样 方法 从 原始 样本 中 提取 多 个 样 


(1) 


表 1 遥感 影像 数据 信息 


Tab.1 Remote sensing image data information 


波段 名 称 Landsat 波 段 ” ”分辨 率 /m 波谱 范围 /nm 波段 名 称 MODIS 波 段 ”分辨 率 /m 波谱 范围 /nm 
Blue 2 30 450~515 Blue 3 500 459~479 
Green 3 30 525~600 Green 4 500 545~565 
Red 4 30 630~680 Red 1 500 620~670 
NIR 5 30 845~885 NIR 2 500 841~876 
SWIRI 6 30 1560~1651 SWIRI 5 500 1628~1652 
SWIR2 7 30 2100~2300 SWIR2 6 500 2105~2155 


i: Blue Green „Red „NIR SSWIR1 SWIR2 分 别 为 蓝光 . 绿 光 、 红 光 . 近 红 外 .中 红外 1 .中 红外 2 波段 。 下 同 。 
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FHEWY 


本 ,为 每 个 自 举 样本 构建 决策 树 ,然后 将 多 个 决策 
树 的 预测 合并 为 最 终 预测 结果 的 平均 值 , 它 构造 的 
不 同 训练 集 可 以 增加 回归 模型 之 间 的 差异 ,从 而 提 
高 组 合 回归 模型 的 外 推 预测 能 力 ,因此 在 土壤 盐 滇 
化 预测 方面 具有 广泛 应 用 ,并 且 效 果 良 好 ,本 文通 
过 Python 3.7 实 现 该 模型 的 算法 过 程 。 本 研究 按照 
7:3 对 样本 进行 划分 ,70% 作 为 预测 集 ,30% 作 为 验 
证 集 。 在 模型 构建 中 ,各 个 特征 变量 对 EC 值 的 贡 
献 率 存在 差异 '”。 本 文 基于 影像 数据 提取 了 6 个 特 
征 波段 反射 率 、 以 及 应 用 范围 较 广 的 植被 指数 和 盐 
分 指数 ,以 EC 值 作为 目标 值 ,应 用 极端 梯度 提升 5 
算法 对 提取 的 指数 进行 特征 重要 性 筛选 ,优选 出 重 
要 性 高 于 0.01 的 6 个 植被 指数 和 9 个 盐分 指数 参与 
建 模 ( 表 2)。 


3 结果 与 分 析 


3.1 融合 结果 与 精度 评估 

对 融合 影像 和 验证 影像 进行 假 彩色 合成 ,以 假 
彩色 显示 ,在 目 视 解 译 的 基础 上 ,对 比分 析 融 合影 
像 和 验证 影像 可 知 (图 2) , ESTARFM 模型 FSDAF 
模型 的 时 空 融合 效果 良好 ,各 种 地 物 大 体 上 纹理 清 
晰 ,图 中 所 显示 的 水 体 农田 .道路 轮廓 特征 清晰 准 
确 ,与 原始 影像 在 空间 上 基本 保持 一 致 。 为 了 清楚 
观察 融合 影像 空间 细节 ,将 融合 影像 的 放大 到 像素 
级 ,对 比 2 幅 融 合影 像 和 验证 影像 细节 差异 可 知 ， 


FSDAF 融合 影像 对 水 体 反映 效果 较 明 显 。 

为 进一步 验证 两 种 融合 模型 的 精度 ,分 别 对 两 
种 模型 融合 影像 的 6 个 特征 波段 反射 率 与 Landsat 
验证 影像 对 应 波段 反射 率 进行 对 比分 析 ( 表 3) ,并 
选取 融合 一 致 性 程度 较 高 的 2 个 波段 (Red .SWIR2) 
进行 2D 散 点 图 显示 (图 3)。 

从 表 3 可 以 看 出 ,ESTARFM 模 型 的 融合 精度 展 
高 于 模型 FSDAF; 两 种 融合 模型 均 是 SWIR2 波 段 融 
合 效 果 最 好 ,SWIR1 波段 融合 效果 最 差 ,融合 结 
SRA) NT AAR , ESTARFM FAHY R°?(SWIR2 )=0.8444, 
R (SWIR1) =0.6358; FSDAF 模型 的 尼 (SWIR2) = 
0.7493 , R° (SWIR1 )=0.4909, ESTARFM 融合 影像 的 
各 个 特征 波段 一 致 性 程度 均 高 于 FSDAF 融合 影像 ， 
其 中 ESTARFM 融合 影像 的 Red .Green Blue Y Et mi 
合 精度 尼 分 别 为 0.8066 .0.6996 .0.6443;FSDAF 融合 
影像 的 Red .Green Blue 波段 融合 精度 尺 分 别 为 
0.6999 ,0.5525 .0.4959。 

ESTARFM 模型 的 融合 影像 中 ,对 土壤 盐 渍 化 

言 息 较 为 敏感 的 Red 波段 .SWIR2 波段 反射 率 与 
Landsat 验证 影像 对 应 的 波段 反射 率 的 一 致 性 尼 分 
别 达 到 0.8066 .0.8444( 图 3) ,2D 散 点 图 结果 显示 和 良 
好 ;在 FSDAF 模 型 的 融合 影像 中 ,对 土壤 盐 渍 化 信 
息 较为 敏感 的 Red 波 段 .SWIR2 波 段 反射 率 与 Land- 
sat 验证 影像 对 应 的 波段 反射 率 的 一 致 性 尼 分 别 达 
到 0.6999 .0.7493( 图 4) ,融合 效果 较 次 于 ESTARFM 


表 2 参与 反 演 的 变量 指数 及 计算 公式 


Tab.2 Variable indices and calculation formulas involved in inversion 


指数 


公式 


归 一 化 植被 指数 (NDVD" 
增强 型 植被 指数 (EVD 

扩展 差 值 植被 指数 (EDVD"" 

比值 植被 指数 (RVD™ 

扩展 的 归 一 化 植被 指数 (ENDVD"™ 
扩展 的 增强 型 植被 指数 (EEVD™ 
盐分 指数 (SIT)P 

SD™ 

sn)" 

S12)" 

S13)" 

归 一 化 盐分 指数 (NDSDP 

盐分 FE oy $1) 
盐分 指数 (S2)P9 
盐分 指数 (S3)P9 


Bt 
> 
mcs mes moe me m 


Fe ee 


pif 


pail 


NIR-Red)/(NIR+Red) 
2.5x{(NIR-Red)/(NIR+6Red-7.5Blue+1)] 
IR+SWIR1-Red 

NIR/Green 
NIR+SWIR2—Red)/(NIR+SWIR2+R ed) 
2.5x[(NIR+SWIR 1—Red)/(NIR+SWIR 1+6Red-7.5Blue+1)] 
Red/NIR)x100 

BluexRed)”* 

GreenxRed)”” 

Green*+Red*+NIR’)”* 

Red+Green 

Red-NIR)/(Red+NIR) 

Blue/Red 

Blue—Red)/(Blue+Red) 

GreenxRed)/Blue 
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(a) ESTARFM 融 合影 像 (d) ESTARFM 融 合影 像 细 节 图 


0 5 km 


L 


(b) FSDAF 融 合影 像 (e) FSDAF 融 合影 像 细 节 图 


0 5 km 


[wm | 


(c) Landsat 影 像 , ( Landsat 影 像 细节 图 


图 2 融合 影像 与 验证 影像 合成 假 彩色 局 部 对 比 图 


Fig. 2 Partial comparison of false color synthesized from fusion image and verification image 


RI 融合 影像 与 验证 影像 波段 反射 率 特征 分 析 


Tab.3 Analysis of the reflectance characteristics of the fusion image and the verification image band 


模型 Blue Green Red NIR SWIR1 SWIR2 
ESTARFM — y=-128.8+1.0848x  y=-113.7+1.0552x  y=-272.8+1.1213x  y=-91.10+1.0880x  y=407.0740.8961x  y=-59.69+1.0485x 
模型 民 =0.6443 R’=0.6996 R’=0.8066 R’=0.7496 R’=0.6358 R’=0.8444 
FSDAF y=-263.2+1.2043x y=-411.1+1.2268x y=-521.9+1.2656x y=-11.97+1.1058x y=445.28+0.9519x y=-267.6+1.1818x 
模型 R’=0.4959 R’=0.5525 R’=0.6999 R’=0.6202 R’=0.4909 R’=0.7493 


注 :* 为 验证 影像 波段 反射 率 ;y 为 融合 影像 波段 反射 率 ; 尼 为 融合 精度 。 


3000 「 Ca) Red 波段 5000 F Cb) SWIR2 波 段 
w y=—59.69+1.0485x 
W ”6000 上  y=-272.8+1.1213x Re 4000 H 7 R20.8444 
a R°=0.8066 D 
anh 3000 
B 4000 È 
2 = 
& & 2000 
: : 
I 
是 2000 
1000 1500 2000 2500 3000 3500 4000 1000 2000 3000 4000 5000 
验证 影像 Red 波 段 反 射 率 验证 影像 SWIR2 波 段 反 射 率 


注 :Red 为 红 光 波段 ;SWIR2 为 中 红外 2 波段 。 下 同 。 
图 3 ESTARFM 模 型 融合 影像 与 验证 影像 2D 散 点 图 
Fig. 3 2D scatter plot of ESTARFM model fusion image and verification image 
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8000 r (a) Red 波段 
y=—521.9+1.2656x 
6000 R?=0.6999 


4000 


2000 


融合 影像 Red 波 段 反射 率 


VAY L 5 L 1 L 1 
1000 1500 2000 2500 3000 3500 4000 
验证 影像 Red 波 段 反射 率 


图 4 FSDAF 模型 融合 影 


像 与 验证 影像 2D 散 点 图 


5000 F (b) SWIR2 波 段 


y=-267.6+1.1818x 
R?=0.7493 


融合 影像 SWIR2 波 段 反 射 率 


1000 2000 3000 4000 5000 
验证 影像 SWIR2 波 段 反射 率 


Fig. 4 2D scatter plot of FSDAF model fusion image and verification image 


模型 的 融合 效果 。 综 合 分 析 可 知 ,ESTARFM 模型 
的 6 个 特征 波段 的 融合 结果 精度 尼 均 在 0.6 以 上 , 相 
较 于 FSDAF 模 型 的 融合 效果 更 优 。 
3.2 不 同 EC 值 的 光谱 反射 率 

为 了 可 视 化 实测 土壤 EC 值 与 时 空 融合 影像 各 
波段 反射 率 之 间 的 关系 ,根据 土壤 盐 涡 化 程度 划分 
标准 ,绘制 出 不 同 盐 度 水 平 区 间 与 对 应 光谱 反射 率 
均值 的 光谱 曲线 关系 图 (图 5)。 对 比 两 种 时 空 融合 
影像 在 不 同 盐 度 的 光谱 曲线 图 (图 5) 可 知 , 针 对 ES- 
TARFM 融合 影像 , 随 着 EC 值 升 高 ,土壤 光谱 反射 率 
也 在 升 高 ,各 波段 光谱 反射 率 曲 线 变 化 趋势 一 致 ， 
从 可 见 光 Red 波 段 开始 陡然 上 升 ,在 NIR 波段 处 达 
到 一 个 峰值 ; 当 EC 值 在 0~2 mS. em! (JERR EA 
壤 ) 之 间 时 ,光谱 反射 率 最 低 , 当 EC 值 在 2~4 mS- cm” 
( 轻 度 盐 汪 化) 之 间 时 ,光谱 曲线 从 Red 波段 开 始 陡 
然 升 高 ,并 且 从 可 见 光 Red 波 段 到 SWIR2 的 反射 率 
一 直 处 于 最 高 状态 。 


FSDAF 融合 影像 , 当 EC 值 在 0~2 mS- cem (GEER 
渍 化 土壤 ) 之 间 时 ,光谱 反射 率 最 低 ,省 EC 值 >16 
ms'cm:( 盐 十 ) 时 ,光谱 反射 率 最 高 ,并且 光 谱 曲 线 
在 SWIR1 波 段 处 达到 峰值 ;其 余波 段 的 光谱 曲线 从 
可 见 光 Blue 波 段 开始 ,反射 率 不 断 升 高 ,在 NIR 波 
段 处 ,反射 率 达 到 峰值 ,这 与 ESTARFM 时空 融合 影 
像 的 变化 趋势 一 致 。 此 外 ,ESTARFM 时 空 融合 影 
像 和 FSDAF 时 空 融 合影 像 中 的 NIR 波 段 处 ,不 同 
EC 值 水 平 的 光谱 特征 不 明显 ,这 一 点 与 前 人 研究 结 
果 一 致 。 
3.3 预测 模型 精度 评估 
土壤 电导 率 (EC) 能 够 有 效 地 反映 土壤 盐 渍 化 
程度 ,在 整个 研究 区 实测 样本 中 ,EC 值 变化 范围 较 
大 , 介 于 0.0476~108.7 mS.ecm! 之 间 , 建 模 输入 样本 
数据 中 ,剔除 2 个 异常 值 , 样 本 70% 作 为 预测 集 (z= 
46) ,30% 作 为 验证 集 (n=21)。 本 研究 共有 6 个 植被 
BAX .9 个 盐分 指数 .6 个 特征 波段 参与 建 模 ,采用 10 


EC/mS:cm! -@ <2 -2~4 + 48 > 816 六 216 


0.5 F (a) ESTARFM 融合 影像 
0.4 + 
= 02 | r 
R 六 
0.1 F 
0.0 1 1 1 1 1 J 


Blue Green Red NIR SWIR1 SWIR2 
波段 


05 - O) FSDAF 融 合影 像 
0.4 H 
0.3 + 
0.2 + 
0.1 上 


0.0 1 1 1 1 1 J 
Blue Green Red NIR SWIRI SWIR2 


波段 


注 :Blue .Green 分 别 为 蓝光 AO CUE BE ; NIR 为 近 红 外 波段 ;SWIR1 为 中 红外 1 波段 。 


图 5 不 同 盐 度 土壤 在 时 空 


融合 影像 中 反射 光谱 曲线 


Fig. 5 Reflectance spectrum curves of soils with different salinities in spatio-temporal fusion images 
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h 120 F(a) ESTARFM 了 融合 影像 T 1297 (o) FSDAF 融 合影 像 
© 100 5 100 
gq | 970.5405x+7.6497 A g y=0.5098x+7.7588 
£ R?=0.9268 & R°=0.8987 
@ 60 = 
= = 
5 40 E 
3 20 Q 
0 
7 120 F (c) Landsat 影 像 
= 100 
80 y=0.5662x+6.5255 


R=0.9103 3 
®. 


EC 实测 值 /mS:cmrl 


图 6 基于 3 种 影像 EC 值 拟 合 
Fig.6 EC values fitting based on 3 kinds of image 


倍 交叉 验证 方法 进行 模型 验证 。 为 评估 预测 模型 
的 性 能 ,选用 尼 作 为 预测 建 模 的 评价 指标 , 尼 值 越 
大 ,模型 精度 越 高 。 

对 比分 析 融 合影 像 验证 影像 的 EC 值 建 模 , 结 
果 显 示 ( 图 6) ,基于 时 空 融合 影像 构建 的 EC 预测 模 
型 精度 良好 ,其 中 ,基于 ESTARFM 融合 影像 构建 的 
EC 值 预 测 模型 效果 最 好 (六 =0.9268 ) ; 基于 Landsat 
影像 次 之 (R=0.9103) ;FSDAF 融合 影像 构建 的 EC 
值 预测 模型 精度 最 低 ( 尼 =0.8987)。 预 测 模型 结果 
理想 ,说 明 选 用 的 两 种 时 空 融合 模型 对 MODIS 和 
Landsat 影像 进行 时 空 融合 ,能 够 获取 高 时 空 分 辩 率 
影像 ,并 在 EC 值 预测 模型 构建 中 应 用 效果 良好 。 


4 讨论 


地 表 与 大 气 环境 需要 实时 精细 的 遥感 监测 ,高 
时 空 分 辩 率 的 遥感 影像 又 较 难 获取 ,时 空 遥 感 影像 
融合 技术 具有 成 本 低 .可 靠 性 强 等 优点 ,近年 来 取 
得 了 较 大 的 研究 进展 ,在 植被 长 势 " 洪 水 监测 中 、 
温度 监测 "中 等 方面 已 取得 一 定 的 研究 成 果 。 本 研 
究 选 用 两 种 较为 成 熟 的 时 空 融 合算 法 ,对 Landsat 和 
MODIS 遥感 影像 进行 时 空 融 合 ,使 融合 影像 兼 具 
Landsat 数 据 的 空间 特征 和 MODIS 数据 的 时 间 特 征 ， 
而 现今 时 空 融 合算 法 都 无 法 准确 预测 土地 覆盖 类 
型 的 微小 变化 ,因此 时 空 融 合 后 的 影像 在 空间 特征 
上 与 Landsat 原 始 影像 产生 差异 ,只 能 不 断 改进 时 空 


融合 算法 来 减 小 差异 。 而 ESTARFM 模型 是 在 
STARFM 模型 基础 上 进行 改进 ,提升 了 对 异 质 性 区 
域 预 测 的 精度 和 适用 性 ,人 研究 显示 ,ESTARFM 模型 
对 盐 涡 化 特征 波段 的 时 空 融合 精度 优 于 FSDAF 模 
型 ,对 土壤 盐分 的 拟 合 精度 也 更 高 。FSDAF 模型 结 
合 混合 像 元 分 解 ,加 入 薄板 样 条 函数 捅 值 .距离 加 
A ,能 准确 预测 土地 覆盖 类 型 发 生 突变 的 区 域 ,已 
有 研究 显示 FSDAF 算 法 对 洪水 预测 效果 良好 "” ,本 
人 研究 中 土壤 盐 涡 化 是 缓慢 发 生 过 程 ,对 FSDAF 算 法 
的 优势 体现 不 明显 ,因此 基于 ESTARFM 融合 影像 
的 土壤 盐分 拟 合 效果 更 好 ,ESTARFM 融合 影像 对 
土壤 盐 渍 化 预测 更 具 适 用 性 。 此 外 ,土壤 盐 渍 化 发 
生 过 程 复 杂 ,为 确保 采样 获取 的 “点 状 " 信 息 能 代表 
土壤 盐 涡 化 分 布 的 “ 面 状 ”信息 ,本 研究 进行 设计 采 
样 时 ,借助 MODIS 数 据 , 以 及 DEM 衍 生变 量 (坡度 
和 地 形 湿 度 指数 ) 和 土壤 质地 数据 ,基于 代表 性 等 
级 设计 采样 ,每 个 采样 点 样 方 为 30 mx30 m, 使 用 五 
点 采样 法 采集 土壤 样品 ,并 选取 与 采样 时 间 相 对 应 
的 遥感 影像 , 重 采样 为 30 m 分 辨 率 后 再 进行 时 空 融 
合 ,使 采样 点 的 样 方 范围 与 遥感 影像 的 像素 相 匹 
配 ,确保 土壤 盐 涡 化 的 预测 精度 。 在 基于 融合 影像 
预测 的 非 盐 溃 土 EC 均值 为 0.38 ms'cm ,面积 占 比 
为 8.96%; 轻 度 盐 渍 土 EC 均值 为 2.61 mS:cm ,面积 
占 比 为 37.20% ; 中 度 盐 渍 土 EC 值 为 5.89 mS .cm ， 
面积 占 比 为 12.09% ; 重度 盐 涡 土 EC 均值 为 11.83 
ms.enm: ,面积 占 比 为 9.55% ; 盐 士 EC 均值 为 31.19 
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mS- cm~, HAR h LEW 32.20%. in EA RI E WF 
究 区 土壤 EC 值 预测 结果 显示 , 非 盐 渍 土 的 EC 均值 
为 0.41 mS .cm , 轻 度 盐 渍 化 土壤 EC 均值 为 2.45 
ms.cem' ,中 度 盐 渍 化 土壤 EC 均值 为 6.25 mS-cm"', 
重度 盐 渍 化 土壤 EC 均值 为 12.07 mS- cm”, $R EEC 
均值 为 40.25 mS .cm … ,经 非 参数 检验 ,其 结果 与 本 
文 预测 结果 并 没有 显著 差异 。 

时 空 融合 影像 在 土壤 盐 渍 化 研究 中 具有 很 大 
应 用 空间 ,Han 等 ”基于 时 空 融合 影像 探究 突 发 降 
水 事件 对 土壤 盐分 预测 的 影响 。 本 文 研究 利用 目 
前 较为 成 熟 的 时 空 融 合算 法 应 用 于 土壤 盐 渍 化 预 
测 的 适用 性 ,但 不 同时 空 融合 算法 的 优势 有 待 控 
掘 , 需 要 进一步 寻求 对 土壤 盐 渍 化 特征 提取 更 具 优 
势 的 时 空 融合 算法 ,提升 遥感 影像 在 土壤 盐 渍 化 预 
测 中 的 应 用 价值 。 


5 结论 

(1) 应 用 ESTARFM 模 型 TSDAF 模 型 能 够 较 好 
的 融合 出 基于 MODIS 和 Landsat 遥感 影像 的 高 时 空 
分 辩 率 影像 ,从 目 视 解 译 结果 来 看 ,融合 影像 空间 
细节 信息 比较 清晰 ,与 验证 影像 空间 信息 基本 一 
致 ,说 明 融 合 结果 理想 。 不 同 模型 对 于 影像 各 个 波 
段 的 融合 结果 均 有 差异 ,这 可 能 是 由 不 同 模型 的 精 
度 误差 所 造成 1。 

(2) 两 种 时 空 融 合影 像 在 不 同 EC 值 水 平 的 光 
谱 特 征 曲 线 变 化 明显 ,基于 ESTARFM 融合 影像 构 
建 的 EC 值 预测 模型 精度 最 好 , 记 =0.9268; 基 于 FS- 
DAF 融合 影像 构建 的 EC 值 预测 模型 精度 良好 ,R= 
0.8987 ; 基于 Landsat 验证 影像 构建 的 EC 值 预 测 模 
型 精度 较 好 ,R=0.9103。 由 此 说 明 时 空 融 合 技术 在 
土壤 盐 渍 化 遥感 监测 中 具有 很 好 的 应 用 价值 和 发 
展 空间 ,对 比 两 种 常用 的 时 空 融合 影像 在 盐分 预测 
方面 的 应 用 效果 ,基于 ESTARFM 的 融合 影像 更 适 
用 于 盐分 预测 。 不 同 环境 协 变量 对 土壤 盐 溃 化 遥 
感 监测 精度 都 有 影响 ,已 有 学 者 尝试 选用 不 同 的 
建 模 因子 .不 同 的 模型 进行 土壤 盐分 预测 ” 。 

(3) 基于 MODIS 和 Landsat 的 时 空 融合 影像 能 
有 效 提升 遥感 影像 的 时 空 分 辩 率 ,解决 时 相 数据 缺 
失 问 题 , 对 比 ESTARFM 模型 和 FSDAF 模 型 的 时 空 
融合 结果 可 知 ,ESTARFM 模型 能 够 融合 出 精度 更 
高 的 时 空 分 辨 率 影像 。 在 土壤 盐 渍 化 遥感 监测 中 ， 


深入 挖掘 时 空 融 合影 像 在 土壤 盐 渍 化 遥感 监测 方 
面 的 应 用 前 景 ,有 助 于 提升 遥感 数据 的 利用 价值 ， 
提升 盐 渍 化 遥感 监测 在 实际 应 用 中 的 经 济 效益 和 
实用 性 。 另 外 ,本 文选 用 GEE 平 台 获 取 影 像 数据 ， 
平台 在 线 完成 预 处 理 过 程 ,高 效 快捷 ,提升 数据 预 
人 处理 速度 ,也 说 明 利用 谷歌 云 计算 平台 获取 逐 感 影 
像 ,进行 土壤 盐 溃 化 监测 同样 可 以 提升 效率 ,为 土 
坏 盐 渍 化 遥感 监测 提供 了 与 时 俱 进 的 利用 价值 和 
发 展 空间 。 
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Exploring the application of MODIS and Landsat spatiotemporal fusion images 
in soil salinization: A case of Ugan River-Kuqa River Delta Oasis 


ZHAO Qiaozhen'’, DING Jianli™, HAN Lijing”, JIN Xiaoye'’, HAO Jianping‘ 

(1. Key Laboratory of Smart City and Environmental Modeling of Autonomous Region Universities, School of Resources and 
Environmental Sciences, Xinjiang University, Urumqi 830046, Xinjiang, China; 2. Key Laboratory of Oasis Ecology, Ministry of 
Education, Xinjiang University, Urumqi 830046, Xinjiang, China; 3. Ministry of Natural Resources, Central Asia Geographic 
Information Technology Development and Application Engineering Technology Innovation Center, Urumqi 830046, Xinjiang, 
China; 4. Ugan River Basin Management Office, Aksu Prefecture, Xinjiang, Aksu 842000, Xinjiang, China) 


Abstract: Remote sensing monitoring of soil salinization relies on high spatiotemporal resolution images, which 
are difficult to obtain due to high cost, revisit period of the satellite, and weather effects. The spatiotemporal 
fusion technology of remote sensing images has the advantages of low cost and good reliability and has been 
widely used in several studies, such as vegetation growth, flood monitoring, and temperature monitoring. This 
study investigates the effectiveness of spatiotemporal fusion images in soil salinization monitoring at Ugan 
River- Kuga River Delta Oasis, south Xinjiang, China. An enhanced spatial and temporal adaptive reflectance 
fusion model and flexible spatiotemporal data fusion model were employed to fuse MODIS and Landsat images 
and generate high spatiotemporal resolution images. A random forest prediction model for soil conductivity (EC) 
was created using the fusion images. Furthermore, the precision of both models was evaluated. Results show that 
the fusion images made by the ESTARFM and FSDAF models have high spatial resolution and good uniformity 
compared with the verification image. Verification revealed the reliability of the fusion results from the 
ESTARFM algorithm, R’?(Red)=0.8066 and R’(SWIR2)=0.8444, and the FSDAF algorithm, R’(Red)=0.6999 and 
R*(SWIR2)=0.7493. The EC prediction model created using the ESTARFM fusion images had the best accuracy 
with R° of 0.9268. Meanwhile, the R? of the EC prediction model constructed from FSDAF was 0.8987 and that 
constructed based on Landsat verification images was 0.9103. This finding showed that spatiotemporal fusion 
technology has good application in the remote sensing monitoring of soil salinization. Comparison of the salinity 
prediction results from the two spatiotemporal image fusion models revealed that the ESTARFM model is 
suitable for salinity prediction. 


Key words: temporal and spatial fusion; random forest; electric conductivity; soil salinization 


